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ABSTRACT 

Near-infrared H2 1-0 S(l) and Br 7 velocity curves along the major 
axis of NGC 253 have revealed a central velocity gradient that is seven 
times steeper than that shown by the optical velocity curve. This 
is interpreted as an optical depth effect due to dust. Approximately 
35 mag of visual extinction in the center is required to match the SW 
side of the optical velocity curve. The spatial variation of the ratio of 
these lines to the 12 CO ( J = 1 — >• 0) line is compared among starburst 
galaxies NGC 253, M82, and NGC 4945 to investigate the excitation 
mechanism responsible for the H2 1-0 S(l) line. 

Subject headings: galaxies: individual (NGC 253) — galaxies: kine- 
matics and dynamics — galaxies: starburst — infrared: galaxies — 
ISM: dust, extinction 
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1. INTRODUCTION 

The recent development of near-infrared 
arrays has made it possible to observe the 
kinematics of spiral galaxies using near-infrared 
emission lines. This is important for measur- 
ing kinematics of the nuclear regions m galax- 
ies with heavy dust obscuration since dust, 
which is concentrated in the galactic plane, 
extinguished the visible light emitted near the 
center. At visible wavelengths, only the line 
of sight velocity of the outer regions of the 
galaxy can be seen. There yet few pa- 

pers in the literature with near-IR kinematic 
studies of galaxies. McKeith et al. (1993) es- 
tablished the wavelength dependence of ob- 
served kinematics for the starburst galaxy 
M82, where they found a much steeper veloc- 
ity gradient for the interstellar [S III] A 9069 A 
and stellar Call IR triplet than for the visible 
and near-UV lines. This effect is attributable 
to dust obscuration, as observed also in the 
dusty inclined galaxy NGC 2146 (Prada et 
al. 1994). 

In this paper, we report near-IR observa- 
tions of the nearby starburst spiral galaxy 
NGC 253. The galaxy has already been ob- 
served along the major axis in optical emis- 
sion lines (Ulrich 1978; Munoz-Tunon et al. 
1993), and in the J = 1 — > transition of 
molecular CO (Canzian et al. 1988). There 
is a significant difference between these two 
data sets. To confirm that this difference is 
caused by the presence of dust, we have ob- 
served NGC 253 in the near-IR, where extinc- 
tion is approximately one order lower than in 
the visible spectral region (Rieke & Lebofsky 
1985). 



2. OBSERVATIONS AND DATA RE- 
DUCTION 

Long slit spectra of NGC 253 were ob- 
tained with the IRSPEC spectrograph (Gre- 
del & Moorwood 1991) at the ESO New Tech- 
nology Telescope at La Silla in March 1994. 
A 58 X 68 SBRC InSb array was used as de- 
tector with a pixel size of 2'/2 . The slit was 
120" long and 4'/4 wide. Seeing was about 
1"; convolved with the pixel size this gives 
a FWHM of 2V8 for a point source. A dis- 
tance of 3.4 Mpc (Sandage & Tammann 1975) 
for NGC 253 yields a scale of 16.5 pc arcsec -1 . 
Two slit positions, one on the major axis (p. a. 
52°) and one on the minor axis (p. a. —38°) 
were observed separately at Br 7 A2.1661 fim 
and H 2 1-0 S(l) A2.1218//m. The spectral 
dispersion was 65 km s -1 pix -1 with a veloc- 
ity resolution of 140 km s -1 . Individual obser- 
vations of one minute were interleaved with 
sky frames in the follow sequence: object-sky- 
sky-object. In total, 14 minutes integrations 
were made on the major axis and 8 minutes 
on the minor, in each spectral range. Fig- 
ures 1 and 2 show the Br 7 and H 2 intensity 
profiles at different positions along the major 
axis. The intensity profiles of both species 
along the minor axis are also plotted in Fig- 
ures 3 and 4. In table 1,2,3,4 are listed the 
Br7 and H 2 line fluxes and LSR radial ve- 
locities at different positions taken from the 
major and minor axis spectra of NCG 253. 

Flatfield and spatial distortion correction, 
sky subtraction, and wavelength and flux cal- 
ibration were performed using the IRSPEC 
package within MIDAS, developed by ESO. 
Flux calibration was achieved by interpolat- 
ing photometric data from the standard stars 
HR5056 and HR2294, in the absence of spec- 
trophotometric data. 

The radial velocities were measured by fit- 
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ting Gaussian profiles using the IRAF task 
FITLINES developed by J. Acosta-Pulido at 
the I AC. Emission line fluxes were determined 
by adding pixels. Velocity measurement er- 
rors combine photon noise, the uncertainty in 
the continuum level and the overall goodness 
of fit. The velocity measurements were cor- 
rected to the local standard of rest (LSR). 

3. RESULTS 

3.1. Near-IR circumnuclear kinemat- 
ics 

We derived Br 7 A2.1661 /mi and H 2 A2.f228 /mi 
position- velocity curves along the major axis 
of the highly obscured nuclear region of NGC 253 
(Figure 5a). The zero position for the ve- 
locity curve is chosen from the continuum 
peak. We measure a systemic velocity Vlsr = 
235 ± 11 kms -1 , in agreement with Vlsr = 
229 km s -1 from CO data (Canzian et al. 
1988) and the optical velocity 248 ± 10 km s" 1 
(Ulrich 1978). On the SW side of the nu- 
cleus, the two near-IR rotation curves are 
seven times steeper than the [N II] A6584A 
optical rotation curve from Ulrich (1978) (see 
Figure 5a). We attribute this difference to 
heavy dust extinction to the SW of the nu- 
cleus (Rieke & Low 1975; Scoville et al. 1985; 
Walker et al. 1988; Pina et al. 1992; Sams et 
al. 1994). A [Nil] and Ha position-velocity 
curve with better spatial and velocity resolu- 
tion (Munoz-Tunon et al. 1993) shows several 
dips indicative of extinction near the nucleus. 

The 12 CO (J = 1 -> 0) 115 GHz position- 
velocity curve (Canzian et al. 1988) is su- 
perposed on our near-IR data in Figure 5b 
for the inner 40" . It agrees well with the 
near-IR kinematics, confirming that the wave- 
length dependence of the kinematics is due 
to dust extinction. We estimate from the 



0.6kms _1 pc _1 slope in Br 7, assuming cir- 
cular motions, a dynamical mass within the 
central 20" diameter region of 4 X 10 9 M & (50 
times greater than would have been obtained 
using the [Nil] lines). 

The velocity dispersions at the nucleus po- 
sition are 63± 4 km s -1 and 61± 4 km s _1 , 
for Br7 and H 2 respectively (corrected for in- 
strumental broadening). This is comparable 
to the velocity dispersion of the gas in the nu- 
cleus of M82, determined from the Br7 2.16 
/mi emission line (Gaffney et al. 1995). 

Along the minor axis of NGC 253, it has 
been observed that optical lines are split- 
ted, with the velocity separation ranging from 
about 250 km s -1 at position nearest the nu- 
cleus to about 450 km s -1 further out. This 
splitting has been attributed to shock excita- 
tion by a starburst-driven superwind at the 
surface of a cone or bubble by Heckman et 
al. (1990) and Schulz & Wegner (1992). In 
the inner ±10 " there are no indications of ve- 
locity splitting in optical lines, in agreement 
with the fact that no splitting is seen in the 
infrared Br7 and H 2 intensity profiles along 
the minor axis shown in Figures 3 and 4. 

3.2. Modeling the dusty rotation curves 

A quantitative numerical model is applied 
in order to investigate how dust could be re- 
sponsible for the difference between the opti- 
cal and near-IR velocity curves. The model 
is an expanded version ofthe Triplex model 
by Disney et al. (1989). It consider an in- 
clined galaxy with an exponential distribution 
for the light emitting gas, embedded in an ex- 
ponential distribution of dust. 

Along the line of sight (z), the density of 
gas is given by: 

p' g (z) = p g e- h ^e-^ (1) 
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where h 3 , r 3 are the scale height and scale 
lenght of the gas respectively. A similar ex- 
pression is assumed for the dust distribution 
where h^, r^ are the scale height and scale 
lenght of the dust respectively. 

At every point, it is assumed that the gas 
nas a velocity distribution which has an aver- 
age velocity given by an input rotation curve, 
a velocity dispersion of 10 km s -1 , and is 
isotropic. The galaxy is supposed to rotate 
cylindrically, and for the intrinsic rotation 
curve we have chosen a curve that rises lin- 
early from the center to 150 km s -1 at 4", and 
remains constant at larger radii. We choose 
the central optical depth in V, and convert 
that to the optical depth at other wavelengths 
using the Galactic Extinction Law (Rieke & 
Lebofsky 1985). For the inclination of the 
galaxy we take an angle of 78.5° (Pence 1981). 
The projected rotation curve is finally deter- 
mined by taking the velocity of the peak in- 
tensity of the line of sight velocity distribution 
at each position along the major axis. 

The scale length of the gas can be taken as 
a free parameter, since it is very difficult to 
measure from the data. Since however very 
few models fit the data, this method provides 
us a way to measure it. The same holds for 
its scale height. We find well-fitting solutions 
for a scale length of 5" and a scale height of 
1". We also have to make a guess for the 
distribution of the dust. A conservative as- 
sumption is that it follows the distribution of 
the stars (e.g. Kaylafis & Bachall 1987) with 
a scale length of 40" and a scale height of 4" 
(Scoville et al. 1985). 

We have plotted the projected rotation 
curve of the model in Figure 5a. We present 
rotation curves corresponding to 0, 2, 20, 30 
and 40 magnitudes of extinction in the galac- 
tic center. The curve for 2 mag of extinction 



fits the near-IR data well, while the [N II] data 
are fitted well, in the SW, by 30-40 mag of 
central extinction in the galactic nucleus. 

Hence, we estimate that approximately 35±5 
mag of visual extinction in the center is re- 
quired to match the SW side of the optical 
rotation curve, which may be compared with 
the value of Ay = 24 mag quoted by Sams et 
al. (1994). We therefore also attribute the 
agreement between the optical and near-IR 
rotation curves in the NE to the presence of 
much less local dust obscuration (Scoville et 
al. 1985). 

3.3. Ionized and Molecular gas distri- 
bution 

The Br 7 emission is strongly peaked at the 
nucleus and we detect it to 12" SW and 18" 
NE along the major axis (Figure 6). The H 2 
1-0 S(l) flux at the nucleus is three times less 
than the Br 7 flux there, while ±8" away from 
the center the H 2 is stronger than the Br 7. 
We do not detect Br 7 further out, while the 
H 2 emission extends to ±20" from the cen- 
ter along the major axis. This distribution 
of the ionized and molecular hydrogen gas in 
NGC 253 is in agreement with Circular Vari- 
able Filter (CVF) spectroscopy with a 19"6 
aperture (Wright & Joseph 1989). We detect 
H 2 emission over an 18" range and Br 7 over 
an 11" range along the minor axis. 

The CO molecular gas emission (Canzian 
et al. 1988) is coexistent with the ionized and 
molecular hydrogen emission along the ma- 
jor axis of NGC 253 (Figure 6). The agree- 
ment of the kinematics of the ionized and ex- 
cited molecular hydrogen with that of the CO 
suggests that the three species coexist in the 
same spatial region. 

The starburst galaxies NGC 253, M82, and 
NGC 4945 have remarkably similar far-IR and 
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radio continuum luminosities (Rieke et al. 
1980; Koornneef 1993; Ghosh et al. 1991). 
All three show two-lobed superwind outflow 
along their minor axes (Heckman et al. 1990; 
Schulz & Wegner 1992) and in all cases the 
H 2 emission is considered to be shock-excited 
by supernova remnants (SNRs) rather than 
by UV fluorescence (for NGC 253 see Mouri 
1994; for M82 see Lester et al. 1990; for 
NGC 4945 see Koornneef 1993 and Moorwood 
& Oliva 1994). 

The spatial variation of the H 2 /Br7 ratio 
however is quite different for the three galax- 
ies (Figure 7). In NGC 253, the ratio is 0.22 in 
the nucleus and grows quickly with increasing 
radius, reaching 2 at 10" radius. While the 
H 2 /Br7 ratio is only about 0.15 in the center 
of M82, it rises to 0.55 at the edge of the star- 
burst core (Lester et al. 1990). In NGC 4945, 
the ratio 2.5 in the central 10" and 

grows to 9 further out (Moorwood & Oliva 
1994). The distribution of the ionized and 
molecular gas in NGC 253 is qualitatively dif- 
ferent from that in the starburst galaxies M82 
(Lester et al. 1990) and NGC 4945 (Koorn- 
neef 1993; Moorwood & Oliva 1994). One 
explanation is that the varied behavior of the 
H 2 /Br7 ratio among the three galaxies is due 
the dependence of H 2 /Br7 on the geometry 
of the emitting regions (Puxley et al. 1990) as 
well as on the evolutionary phase of the star- 
burst (Doyon et al. 1994). Because the ratio 
of H II regions and supernovae is different de- 
pending on the age and geometry of the star- 
burst, a large range in H 2 /Br7 among star- 
burst nuclei is expected (Moorwood & Oliva 
1988, 1990; Puxley et al. 1990; Doyon et al. 
1994). 

In order to study the excitation mechanism 
of H 2 the spatial variation of the intensity ra- 
tio of H 2 A2.1218//m to 12 CO (J = 1 -> 0) 



emission along the major axis of NGC 253 is 
shown in Figure 8. (The H 2 data have been 
convolved to the same resolution as the CO.) 
The ratio decreases sharply with radius, and 
its behavior is quite similar in M82 (CO from 
Lo et al. 1987; H 2 from Lester et al. 1990). 
The ratio has been normalized to unity at the 
nucleus in both cases. Smith et al. (1991) 
have shown the H 2 line strength depends on 
the shock model in molecular clouds. To this 
extent the integrated intensity of the 12 CO 
emission line measures the H 2 column den- 
sity (e.g., Scoville et al. 1987). If the H 2 
is predominently shock-excited then we sug- 
gest that the variation in the H 2 A2.1218/im 
to 12 CO ratio may reflect changed fraction of 
shocked gas due to SNR. If the same excita- 
tion mechanism for the H 2 can be applied to 
NGC 253 and M82, the H 2 /CO ratio (Fig- 
ure 8) shows that shocks may be most com- 
mon near their galactic nuclei. These con- 
clusions however are speculative and model- 
dependent. 

4. CONCLUSIONS 

Near-infrared observations along the ma- 
jor axis of NGC 253 have confirmed that the 
observed optical kinematics are notably af- 
fected by extinction. We attribute the steeper 
central gradient of the infrared velocity curve 
to the smaller line-of-sight optical depth, and 
thus to a deeper penetration into the dusty 
galaxy. Hence, the true mass in the rapidly 
rotating nucleus only becomes apparent at 
infrared or longer wavelengths affected min- 
imally by dust extinction. 

The spatial variations of the H 2 1-0 S(l) 
to 12 CO (J = 1 — > 0) ratio along the ma- 
jor axes of NGC 253 and M82 are similar 
and if the H 2 is considered to be principally 
shocked-excited then the H 2 /CO ratio may 
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trace the fraction of shocked gas due to su- 
pernova remnants, increasing towards their 
nuclei. The spatial variation of the H 2 /Br7 
ratio differs among starburst nuclei, reflect- 
ing possibly differences in the geometries of 
the emitting regions and the ages of the star- 
bursts. 

We thank E.Perez, J.Acosta, L.Cuesta, R.Genzel 
and specially Alicia Gordillo for valuable com- 
ments on the manuscript. F.Prada is the 
holder of a SERC research studentship, at 
QUB, under the IAC-SERC agreement. 
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Figure Captions 

Fig. 1 Br7 intensity profiles at different po- 
sitions along the major axis of NGC 253 

Fig. 2 H 2 2.12 /mi intensity profiles at dif- 
ferent positions along the major axis of NGC 
253 



Fig. 3 Br7 intensity profiles at different po- 
sitions along the minor axis of NGC 253 

Fig. 4 H 2 2.12 /mi intensity profiles at dif- 
ferent positions along the minor axis of NGC 
253 



Fig. 5 (a) Measured radial velocities in Br 7 
A2.1661 /mi and H 2 1-0 S(l) A2.1218 /mi along 
the major axis (p. a. 52°) of NGC 253, com- 
pared with the [N II] A6584 A position- velocity 
curve by Ulrich (1978) and the CO rotation 
curve by Canzian et al. (1988). Model rota- 
tion curves including the effects of dust are 
drawn for a central extinction of 0, 2, 20, 30 
and 40 magnitudes. The steepest curve cor- 
responds to no extinction and the shallowest 
to 40 mag. The portion of the near-IR curve 
with low radial velocity to the NE between 
10" and 25" from the centre, can be attributed 
to a real effect, e.g. the dynamical influence 
of the bar (see Canzian et al. 1988) (b) The 
Br 7 and H 2 position- velocity curves are com- 
pared to the CO rotation curve from Canzian 
et al. (1988). The near-infrared curves show 
a slope of 0.6 km s -1 pc -1 . 1" = 16.5 pc. 
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Fig. 6 Comparison of the distribution of the 
Br 7 and vibrationally excited H 2 A2.1218 fim 
emission with the 12 C0 ( J = 1 — > 0) emission 
(Canzian et al. 1988) along the major axis of 
NGC 253. The CO fluxes are appropriate for 
an aperture of 4'/4 on the minor axis by 2'/8 
on the major axis. 

Fig. 7 The spatial variation of the H 2 /Br7 
ratio along the major axis of NGC 253 com- 
pared to that for two similar galaxies, M82 
(Lester et al. 1990) and NGC 4945 (Moor- 
wood & Oliva 1994). 

Fig. 8 The variation of the H 2 1-0 S(l) to 
12 CO ratio is shown for the major axes of 
NGC 253 and M82. Multiply graphed values 
by 80 for M82 (based on flux in the central 
12") and by 230 for NGC 253 (based on flux 
in the central IR resolution element). M82 is 
3.25 Mpc distant, so the scale is 15.8 pc arcsec - 
for M82 and 16.5 pc arcsec" 1 for NGC 253. 



Table 1: Observed Br7 line flux and LSR ra- 
dial velocity along the major axis of NGC 253 



Position Line Flux Vlsr 
(arcsec) (10~ 14 erg cm -2 s -1 ) (km s -1 ) 

-15.4 0.46±0.02 227.41±8.91 

-13.2 0.59±0.03 206.86±5.45 

-11.0 0.71±0.04 187.83±5.00 

-8.8 1.32±0.07 175.96±5.21 

-6.6 2.65±0.09 175.43±2.27 

-4.4 4.34±0.14 179.42±2.44 

-2.2 6.76±0.13 224.77±2.63 

0.0 12.09±0.31 254.94±2.41 

2.2 13.49±0.34 264.35±3.29 

4.4 5.80±0.19 279.39±4.17 

6.6 1.32±0.03 307.02±7.54 

8.8 0.48±0.06 317.22±12.43 

11.0 0.27±0.04 351.03±19.78 
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Table 2: Observed H 2 line flux and LSR radial 
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of NGC 253 




Table 3: Ob 


served Br7 line flux 


and LSR ra- 










dial velocity along the minor axis 


or in zoo 


Position 


Line Flux 


VlSR 




(arcsec) ( 


10~ 14 erg cm -2 s 


_1 ) (km s -1 


) 








-22.0 


0.33±0.01 


151.13±8. 


.59 


Position 
(arcsec) ( 


Line Flux 
10~ 14 erg cm -2 s -1 ) 


(km s -1 ) 


-19.8 


0.40±0.02 


171.76±8. 


.56 
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.20 


2.2 


2.12±0.12 
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-11.0 
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Table 4: Observed H 2 line flux 
dial velocity along the minor axis 


and LSR ra- 
of NGC 253 
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3.49±0.02 


210.65±5, 
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3.94±0.15 


231.73±5. 
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Position 
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Line Flux 
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.70 
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.26 
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225.23±14.00 
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